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ABSTRACT 

(n: 

Context. Precise and accurate determinations of the atmospheric parameters effective temperature and surface gravity are mandatory 
■ to derive reliable chemical abundances in OB stars. Furthermore, fundamental parameters like distances, masses, radii, luminosities 

can also be derived from the temperature and gravity of the stars. 
^ , Aims. Atmospheric parameters recently determined at high precision with several independent spectroscopic indicators in NLTE, with 

^ ■ typical uncertainties of ~300 K for temperature and of ~0.05 dex for gravity, are employed to calibrate photometric relationships. This 

is in order to investigate whether a faster tool to estimate atmospheric parameters can be provided. 

Methods. Temperatures and gravities of 30 calibrators, i.e. well-studied OB main sequence to giant stars in the solar neighbourhood, 
, are compared to reddening-independent quantities of the Johnson and Stromgren photometric systems, assuming normal reddening. 

In addition, we examine the spectral and luminosity classification of the star sample and compute bolometric corrections. 
Results. Calibrations of temperatures and gravities are proposed for various photometric indices and spectral types. Once the lumi- 
nosity of the stars is well known, effective temperatures can be determined at a precision of ~400 K for luminosity classes III/IV and 
, ~800 K for luminosity class V. Furthermore, surface gravities can reach internal uncertainties as low as ~0.08 dex when using our 

• ■ calibration to the Johnson Q-parameter Similar precision is achieved for gravities derived from the /3-index and the precision is lower 

i-C ■ for both atmospheric parameters when using the Stromgren indices [cl] and [u - b]. In contrast, external uncertainties are larger for 

the Johnson than for the Stromgren calibrations. Our uncertainties are smaller than typical differences among other methods in the 
literature, reaching values up to +2000 K for temperature and ±0.25 dex for gravity, and in extreme cases, +6000 K and ±0.4 dex, re- 
^ ■ spectively. A parameter calibration for sub-spectral types is also proposed. Moreover, we present a new bolometric correction relation 

^ ■ to temperature based on our empirical data, rather than on synthetic grids. 

^ I Conclusions. The photometric calibrations presented here are useful tools to estimate effective temperatures and surface gravities of 

non-supergiant OB stars in a fast manner This is also applicable to some single-line spectroscopic binaries, but caution has to be taken 
for undetected double-lined spectroscopic binaries and single objects with anomalous reddening-law, dubious photometric quantities 
and/or luminosity classes, for which the systematic uncertainties may increase significantly. We recommend to use these calibrations 
^ ' only as a first step of the parameter estimation, with subsequent refinements based on spectroscopy. A larger sample covering more 

00 ■ uniformly the parameter space under consideration will allow refinements to the present calibrations. 



Key words. Stars: fundamental parameters, stars: early-type, stars: atmospheres 



p 

■ 1. Introduction cleosynthesis and chemical evolution, since the CCSN are the 
' primary sources of several chemical elements, e.g. oxygen. As 

^ ; Stars with spectral types between ~09 and B2 belong to the 'nat- they are very luminous (~10''-10^ Lq), this group of stars be- 
ll. . uraF sgectroscopic group 'OB', presenting strong H and Hei came an important subject of astronomical research in the past 
P*" ■ lines ([Jaschek & Jascheck 1987). The lower boundary of this decades 

k>( . group (B2) corresponds to lower mass and upper age limits of ~8 

■ Mq and 30 Myr, respectively, and the occurrence of core collapse 

' supernovae (CCSN). OB stars delineate the spiral structure of New developments allow to obtain better quality observa- 
our Milky Way and play an important role in cosmic evolution, tions and more realistic physical models to analyse data, than 
They dominate the light of entire spiral and irregular galaxies those available in the past decades. Photographic spectrograms 
and influence the interstellar medium, star formation and galac- previously used have been replaced by very high-resolution dig- 
tic evolution through their ultraviolet radiation, stellar winds and ital spectra. Furthermore, larger telescopes allow us to achieve 
supernovae explosions. Furthermore, they are contributors to nu- higher signal-to-noise ratios. And the advent of more com- 
plex computers permit to incorporate a more realistic physical 

~r~I ' ; '. „ J , ^ A ' ■ TT- background into our calculations of synthetic spectra that are 

Based on observations collected at the Centro Astronomico His- , . , . ™ . ^ , 

pano Aleman (CAHA) at Calar Alto, operated jointly by the Max- ^^ed to interpret observations. These improvements have not 

Planck Institut fur Astronomic and the Instituto de Astrofi'sica de only ^ significant eff^ect on the determination of basic proper- 

Andalucfa(CSIC), proposals H2001-2.2-011 and H2005-2.2-016. ties of the stars, such as those in the present work, i.e. effec- 

** Based on observations collected at the European Southern Obser- tive temperatures and surface gravities, but also on other quan- 

vatory, Chile, ESO 074.B-0455(A) and from the ESO Archive. titles that strongly depend on them. Examples are the chemi- 

*** Based on spectral data retrieved from the ELODIE archive at Ob- cal abundances, distances, luminosities, radii, masses, and ages, 

servatoire de Haute-Provence (OHP). Improved determinations of these quantities allow us to better 
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Table 1. Star id, spectral type, photometric quantities , and spectroscopic stellar parameters 
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Johnson Mermilliod ('1991); Morel & Magnenat (1978) for stars #1, 13, 21, 22, 2 3 and 27. S tromgren: 
[Hoffleit & Jaschek ( 1982); Gray & Corballv (2009) for stars #6, 7, 11, 19 and 23; IConti & Alschuler ( 
Luminosity class corrected (see App.lAt. Uncorrected classes are, for star #2; III, #8; III, #10: IV, #19 
Spectral type corrected (see App.lAt. Uncorrected types are, for stars #2, 8, 16, 28 and 29; B2, #5, 10; 
Stars with anomalous reddening law. Ry values are; 5.5 for stars #21 and 22 and 4.2 for star #27. 



Hau ck & Mermilhodl ( fl998l) . 
1971) for star #30. 
IV, #26; V. 

B2.5, #3, 12; BO, #13; B0.5,#22; B0.7. 
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understand the physics of stars (their interior, nucleosynthesis, 
evolution) and the structure and evolution of galaxies. 

In the present work, we mainly focus on the temperature 
and gravity scales of non-supergiant OB stars, i.e. main se- 
quence to giant stars (core H-burning). And we also evaluate 
bolometric corrections in function of temperature. We inves- 
tigate whether our recent improvements made in the acquisi- 
tion of observational data, calculation of more realistic theoret- 
ical spectra and in the spectral analysis methodologies, affect 
the overall effective temperature and surface gravity scales of 
those objects. A quantitative spectroscopic analysis for a rep- 
resentative sample of 30 stars in the solar neighbourhood has 
been performed in Nieva & Przvbillal (l2012| hereafter Paper II 
iNieva & Simon-DfazI (1201 iL Paper ID and iBriquet et all hoi l( 
Paper III), where multiple independent ionization equilibria 
were evaluated iteratively together with fits to Balmer lines in 
NLTE in order to reproduce the whole observed spectrum in 
the optical. The analysis resulted in a unique set of atmospheric 
parameters at high precision and accuracy for each star. Based 
on these stellar parameter determination, we provide here new 
calibrations to determine temperatures and gravities of OB stars 
from photometric quantities and their spectral classification. 

The paper is structured as follows. Section |2] presents the 
star sample and observations, Sect.|3]briefly addresses the tech- 
niques used for stellar parameter determinations in previous pa- 
pers, Sect.|4]discusses the temperature and gravity scales to pho- 
tometric data. In Sect. |5] our photometric scales and spectro- 
scopic parameters are compared to other work from the liter- 
ature. Sect. |6] focus on the spectral type calibration to temper- 
ature and Sect. |7] on the bolometric coiTections as a function 
of temperature. Section |8]presents a summary and conclusions. 
Furthermore, App.|A] addresses the spectral type and luminosity 
class examination and App. |B] the spectroscopic indicators for 
temperature and gravity determinations. 

2. The star sample and observations 

Our star sample covers a slightly broader temperature range than 
the standard definition for OB stars, with spectral types from B3 
to 09. In the context of stellar evolution, the stars are located 
between the Zero Age (ZAMS) and the Terminal Age Main 
Sequence (TAMS). Their temperatures range from ~ 16 000 to 
34 000 K, the surface gravities from ~3.6 to 4.3 in logarithmic 
scale, and the masses from ~5 to 20 Mq. The sample is listed 
in Table [1] where the spectral type, photometric quantities and 
atmospheric parameters derived in the Papers I-III are given. 
Object identifications include HD and HR numbers and names 
of the stars and our own numbering scheme to facilitate easy 
identification in the figures and other tables (i.e. stars analyzed 
in Paper I: #1-20, Paper II: #21-29, Paper III:#30). 

For the whole sample, high-resolution echelle spectra at very 
high signal-to-noise (S IN) ratio - ranging from 250 up to over 
800 in B - and wide wavelength coverage were obtained, either 
by own observations, or from archives. Details on the observa- 
tions and data redu ction can be found in Papers I and III and 
for stars #21-29 in 'Simon -DfazI (1201 0'). Spectra for stars #1- 
6 were obtained with Feros on the ESO 2.2m telescope in La 
Silla/Chile, with a resolving power of R- A/AA « 48 000. Stars 
#7-16 were observed with Foces on the 2.2 m telescope at Calar 
Alto/Spain at a resolution ofR^AO 000. Spectra of stars #17-20 
were extracted from the archive of the Elodie echelle spectro- 
graph mounted on the 1.93 m telescope of the Observatoire de 
Haute-Provence/France at 42 000. Spectra for stars #21-29 
were obtained with Fies mounted at the 2.5 m Nordic Optical 




oq T,. 



Fig. 1. Location of the sample stars in the Tf^g-logg diagram. 
Atmospheric parameters and star identification are according to 
Table [U Evolutionary tracks for non-rotating stars of metallicity 
Z = 0.014 are from .Ekstrom et al., (.2012) . The location of the 
ZAMS is indicated. 



Telescope (NOT) in El Roque de los Muchachos observatory 
on La Palma (Canary Islands, Spain). A spectrum obtained with 
Feros for star #30 was extracted from the ESO Archive. 

We aimed at selecting single stars, as second light from a 
companion of a similar spectral class distorts the ratio of line- 
to continuum-fluxes. Thus, objects in SBl systems with much 
fainter companions (e.g. stars #21 and 26) can therefore still 
qualify as targets, because the contribution of the secondary star 
in optical light (UBV, uvby) is marginal. Individual components 
in a visual binary are also considered as targets (e.g. star #13 ). 
The star sample contains normal, jSCephei and slowly pulsating 
B-starfl Exposures taken during a small fraction of the same 
night were co-added. The analysis will thus yield atmospheric 
parameters like effective temperature and surface gravity valid 
only for the moment of observation. For this sample, we have 
adopted reliable photometric data from the literature, these how- 
ever coiTespond to a different moment of observation than for 
the spectroscopic data, which might be of significance for some 
variable stars. 



3. Stellar parameter determination 

Atmospheric parameters were determined simultaneously via 
line-fitting of synthetic spectra computed in NLTE to the ob- 
served data, as described in Papers I-III. Atmospheric parame- 
ters are the effective temperature Teff , surface gravity log g, mi- 
croturbulence ^, projected rotational velocity v sin / and macro- 
turbulence Details on the synthetic spectra and the quantitative 
analysis are summarized in the following. 



3. 1 . Spectrum synthesis in NLTE 

The analysis was based on a grid of synthetic spectra in NLTE 
computed with the code s Detail and Surface (lGiddingslll98ll : 
iButler & Giddings"1985) on prescribed LTE atmospheric struc- 
tures calculated with Atlas9 (Kurucz 199312). Model atoms for 
H, He, C, N, O, Ne, Mg, Si and Fe for different ionization stages 



' Pulsation periods of ySCephei stars vary between ~3 to 8 hours, 
while those in SPB stars vary between -0.5 to 3 days. 
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were employed. Opacities due to hydrogen and helium and the 
individual metals were considered in NLTE, and line block- 
ing was ac counted for i n LTE via Kurucz' Opacity Distribution 
Functions (lKuruczlll993ah . This approach provides an efficient 
way to compute realistic synthetic spectra when the atmosphere 
is close to LTE. In particular, it is equivalent to full NLTE 
model atmosphere calculations for late O-type to early B-type 
stars with weak winds (.Nieva & Przvbilla.2007; .Przvbilla et aL 
l20TTh . 

3.2. Quantitative spectral analysis 

The spectral analysis follows the methodology introduced by 
LNieva (2007) and applied for H/He and C by Nieva & Przvbilla 
(120071 12008D and for other elements by Przvb illa et al.l (l2008l) . 
Stellar parameters are determined through ionization equilibria 
for different species and Balmer lines fits. Spectroscopic indica- 
tors for effective temperature and surface gravity of the star sam- 
ple are summarized in Table IBTTI Fitted spectra for stars #1,6,7 
and 15 (Paper I) and #9 (Paper II), illustrate examples of global 
fits for different spectral classes. Numerous examples of fits to 
stellar spectral energy distributions (SEDs) can also be found in 
Papers I and III. The reliability of our temperature and gravity 
determinations and the minimisation of systematic uncertainties 
are confirmed by simultaneous agreement of several independent 
indicators: ionization equilibria of different atomic species, fits 
to Balmer lines and fits to the spectral energy distributions from 
the UV to the IR. Furthermore, our surface gravities are also 
confirmed by simultaneous match of ionization equilibria, fits to 
Balmer lines and the agreement of spectroscopic distances with 
Hipparcos distances (see Paper I). Note that every model atom 
for NLTE calculations is independent from each other, hence 
multiple ionization equilibria only render 4-5 independent in- 
dicators for temperature and gravity determinations. 

The location of the (apparently slowly-rotating) sample stars 
in the Teff-log g plane is visualized in Fig. [T] where also evo- 
lutionary tracks for non-rotating stellar models (Ekstromet al. 
l20T2h are shown for comparison. Both, the models and the stars 
have the same metallicity of Z- 0.014, the stars' metallicities be- 
ing determined in Papers I-III. The star sample provides us with 
a good coverage of the main-sequence band, ranging from ~6 to 
20 Mq. 



4. Temperature and gravity scales 

Useful temperature sc ales for early -t ype stars have been pre- 
sented by Napiwotzk i et all (119931) . iDaflon et ai] (119991) and 
iLvubimkov et al.l 020021) , where the temperature can be esti- 
mated from reddening-independent photometric quantities (as- 
suming a valid standard reddening-law). The first work scales 
the temperature to the [u - /jj^-parameter from the narrow-band 
Stromgren M^vy-system. In contrast, the other two employ the 
(^■iiidex from the broad-band Johnson U B V-sy stem. In all three 
cases, the temperature determination is based on other indepen- 
dent ph otometric cahbrations or on fits to spectral energy distri- 
butions, j^ubi^tovetal] (HOO^) proposed, in addition, temper- 
ature scales for two significantly different values of gravity. 

Motivated by the aforementioned studies, we investigate 
whether new temperature and gravity scales to Q, {u - b\, the 



2 lu-b] = (u -b)- l.Sib -y) = ci + 2m, + 0.5 (b - y) 

3 Q = (U - B)-0.72(B-V) 



Stromgren [clf] and /S can be established based on our pre- 
viously derived spectroscopic parameters. Our atmospheric pa- 
rameters derived simultaneously from several independent indi- 
cators and high quality spectral models and observations offer 
a unique set of input data for building such scales: uncertain- 
ties are as low as ~l-2% for effective temperatures and ~10% 
for surface gravities. This is in contrast to uncertainties of ~ 5- 
10% and ~25%, respectively, using standard methods. Typical 
uncertainties for effective temperature given in the literature are 
<3-5%, however these do not take into account the large discrep- 
ancies - up to ~6000 K - that exist between different temperature 
determination methods. Error estimates of up to 40-60% in sur- 
face gravity are also found in the literature (see Sect.|5]for more 
details). Likewise, the most accurate photometric data available 
in the literature are adopted for this study. The compilation of 
photometric data is as homogeneous as possible, however for 
some objects, data are adopted from different sources. In addi- 
tion, re-examinations of spectral types and luminosity classes 
were performed for the sample stars, as described in App. |A] 
Stellar parameters, photometric indices, and their respective un- 
certainties, spectral types and luminosity classes adopted for the 
sample are listed in Table[T] 

In the following, we present empirical relationships of atmo- 
spheric parameters and photometric indices based on data of our 
star sample. Note that the star sample does not cover uniformly 
the parameter space under consideration, therefore the relation- 
ships can be refined when adopting a more complete sample. 
Polynomial fits are chosen to directly compare our results to pre- 
viou s works from[N apiwotzki etaiT([T993l) . lOaflon et alJ 119991) 
and Lvubimk ov et al.l (|2002). A few stars do not follow tightly 
the empirical fits for different reasons that need to be further in- 
vestigated. Moreover, the correlation between effective temper- 
ature and surface gravity could be treated in an iterative manner 
in o rder to find an opt i mal so luti on, hke the method prop osed 
in Moon & Dworetskvl (II984 ') and Na piwotzki et alj (I1993I) . We 
aim at investigating this correlation with a larger sample of stars 
covering uniformely the parameter space in a further paper. 

4.1. Spectroscopic parameters vs. Q 

The reddening-free Johnson Q-index is particularly interesting 
for OB stars, as there are many distant objects in the Galactic spi- 
ral arms that at present do not have measurements in photomet- 
ric systems other than the Johnson-system. Figure |2^ shows the 
relation between our spectroscopic temperatures and Q, while 
Fig-Hh displays the behavior of red / \ogg vs. Q. The data from 
Table [U are indicated by symbols and the polynomial fits by full 
lines, encoded for different luminosity classes as shown in the 
legend (in color in the online version). Four objects of luminos- 
ity class V had to be excluded from the analysis - marked with 
a cross in Fig.|2^ (# 21, 22, 24 and 27). These stars are located 
in the Orion OB association. Stars # 21, 22 and 27 are located in 
a region of dense nebulosity and present anomalous reddening 
(see Table[T]). In these cases, the Q-index should be re-evaluated 
using an appropriate reddening law, which is, however beyond 
the scope of the present paper. 

The distinction of different luminosity classes reveals tighter 
trends of Teff vs. Q than pre viously reported. The tren ds are con- 
sistent with those found bv ' Nieva & Przvbillal (l2008h for a sub- 
set of stars (#1-6), where the parameter determination was based 
only on H, He and C. However, that sample was insufficient for 

* [cl] =cl-0.2(b-y) 

' yS is a measure of the stellar H/3 line intensity 
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Fig. 2. Functional relationships for the dependence of T^a- (left panels) and log g (right panels) on the photometric parameters Q 
(a,b), [u - b] (c,d), [cl] (e) and /3 (f). Symbols correspond to the sample stars and curves are polynomial fits to the data for different 
luminosity classes, encoded according to the legend. A few stars have been excluded from the fits (crosses). Typical error bars are 
indicated. See the text for details. 



a statistical analysis. In the present work, the polynomial fits for 
luminosity classes III and IV converge into a single curve and the 
scale for luminosity class V is ~ 1500 to 2000 K hotter. The min- 
imisation of uncertainties in our parameter determination and the 
high accuracy of the photometric data allow clear trends for dif- 
ferent luminosity classes to be distinguished. In contrast, stan- 
dard parameter determinations may provide larger systematic 
uncertainties (see Sect.|5]l, that would prevent the relations found 
here to be identified. 

The relationship between the temperature and Q for different 
luminosity classes, i.e., different evolutionary stages, displayed 
in Fig. can be expected by the fact that stars closer to the 



ZAMS are more compact and bluer, therefore hotter than the 
more evolved giants at similar masses. 

In contrast, the ratio Teff / log g vs. Q does not show a depen- 
dence on the luminosity class (Fig. Wp), since all three polyno- 
mial fits converge approximately to a single solution. Therefore, 
we adopt one fit to the whole sample (except for the excluded 
stars), indicated by the thick black curve. 

Temperature and gravity scales to the Johnson Q-parameter 
are proposed, based on the curves displayed in Figs. |2^,b. 
Polynomial fits result in the following equations for the differ- 
ent luminosity classes: 

V : Te = (22.59 + 48.57 Q + 70.26 Q^) 10^ K (1) 
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III, IV : = (33.42 + 78.40 Q + 87.19 Q^) 10^ K 



Tq{W, IV or III) 
(6.67 + 13.30 Q + 16.67 Q^) 10^ K 



(2) 



(3) 



The range where these equations are valid are delimited by 
the observational data, as indicated in Figs. |2^,b. Note that the 
parameter range of interest is not uniformly sampled by the ob- 
servations (Fig. [T]), a further addition of few hotter stars with 
lower surface gravities will allow to refine the fits from the 
present work. The internal and external uncertainties of Tq and 
log^Q are listed in Table |2] The internal uncertainty measures 
the Icr-dispersion of the residuals of the fitted values and the 
spectroscopic parameters and it depends on the quality of the 
polynomial fits. The external uncertainty is calculated via error 
propagation of the parameters involved and it depends on the 
uncertainties of the spectroscopic and photometric parameters, 
and on the reliability of luminosity classifications. Luminosity 
classes are the largest sources of external errors , as the relation- 
ship between temperature and the Q-'mdtx for luminosity class 
V is up to ~2000 K hotter than that for III,IV. The precision of 
the parameter determination is as indicated in Table|2]only when 
accurate photometric data and a reliable luminosity class is given 
for the sta:0. 

When the luminosity class is not known or is not reliable, 
the uncertainty in the derived temperature increases to ~1500K, 
which is still smaller than the differences between various tem- 
perature determinations in the literature (they can amount up 
to ~2000K and in extreme cases up to ~6000K; see Sect. |5] 
for more details). In the case of the gravity, the lack of pre- 
cise luminosity classes causes an error of ~0.2 dex, that is also 
smaller than the maximum differences shown in the literature of 
+0.4 dex. Uncertainties in temperature of ±1500K and in grav- 
ity of +0.2 dex are, however, too large for some further appli- 
cations, such as precise chemical abundance determinations or 
precise derivation of mass, radius and distances of stars. Despite 
the luminosity classification can compromise the accuracy of 
temperature and gravity, our calibrations are a useful tool to esti- 
mate such parameters for further applications. Therefore, Tq and 
log gQ can be employed as starting points to be refined using, e.g. 
spectroscopic analyses. 

4.2. Spectroscopic parameters vs. [u - b] 

In this section we investigate a scale based on the narrow-band 
Stromgren colors for estimating effective temperatures and grav- 
ities. Figures |2};,d show the relationship between [u - b] and 
our spectroscopic temperatures and T^ff/ \ogg, respectively. The 
data are encoded like in Figs. |2^,b. Four objects had to be ex- 
cluded from the analysis - marked with a cross in Fig.|2};. Stars 
# 21, 22 and 27 have anomalous reddening and star #30 is appar- 
ently too hot for the relationship found for the rest of the main 
sequence starfl Similarly to the case of the Q-parameter, trends 
of Teft vs. [u - b] are found for different luminosity classes. In 
this case, the polynomial fits for luminosity classes IV and V 
converge into a single relation, while the scale for luminosity 
class III is ~2000 to 4000 K cooler. In addition, Teff/logg vs. 



* Luminosity classes from Sim bad from differ ent sources 

might differ from each o ther. See iNieva & Przvbillal ( l2008h and 
iFim stein & Przvbillal d2012h for corrections of luminosity classes of 
the Simbad database. 

^ iNapiwotzkiet all ( 119931) propose an upper limit of 30 000 K for this 
relationship. 



[u - b] does also not indicate a dependence on the luminosity 
class (Fig.|2}l). The curves displayed in Figs.|2]:,d are employed 
to propose a temperature and gravity scale to the Stromgren pho- 
tometric parameter [u - b]. Polynomial fits to the data result in 
the following equations for the different luminosity classes: 

IV, V : r[„_i,] = (35.82 - 54.86 [u-b]+ 38.07 [u - bf) 10^ K(4) 



III : r[„_i] = (31.45 - 46.00 [u-b] + 32.96 [u - bf) W K (5) 



logg[«-i] = 



r[„_i] (V, IV or III) 



(7.98 - 8.67 [u-b] + 3.48 [u - b]^) 10^ K 



(6) 



Similarly to the case of the Q-parameter, the range where 
these equations are valid are determined by the observational 
data. The internal and external uncertainties of T[,i-h] and 
log g[u-b] are listed in Table |2] When the luminosity class is not 
known or not precise, the error in the derived temperature in- 
creases to ~2000 K and for gravity to ~ 0.3 dex, both larger than 
in the case of the calibrations to Q. 

4.3. Spectroscopic parameters vs. [cl] and /3 

The Stromgren parameters [cl] and /3 have been defined as tools 
to estimate T^ff and \ogg in early-type stars. A single trend for 
the effective temperature vs. [cl] (Fig. |2^) and two curves for 
\ogg vs. y6 (Fig.|2f) are found. For the latter, the relationship for 
luminosity class V is ~0.3dex higher than for III/IV. Stars #4 
and 26 have been excluded from the fits because their gravities 
are too low in comparison with the rest of the stars of luminosity 
class III/IV (they He slightly above the TAMS in Fig.[B. Star # 
21 is also excluded, as there is a large nebular emission feature 
centered on H/3, that causes a significantly reduced value of the 
/3-index. As the parameter range of interest is not uniformly sam- 
pled by the observations (Fig.[TJ, a further addition of few hotter 
stars with low surface gravities and/or cooler stars on the ZAMS 
might produce a flattening of the slope for luminosity class V in 
Fig.|2f, which remains to be investigated. 

Polynomial fits to the present data result in the following 
equations: 



T[ci] = (27.57 - 52.92 [cl] + 61.25 [cl]^) 10^ K 



V: log^ = 8.75 - 1.73yS 



III, IV: \ogg^ 12.14- 3. 14yS 



(7) 



(8) 



(9) 



The uncertainties are listed in Table|2] If the luminosity class 
is not well known, the error in the surface gravity determination 
increases to ~0.25 dex. 



4.4. Spectroscopic parameters vs. fitted values 

Internal uncertainties of our calibrations are analysed here via 
residuals of the fitted values to the spectroscopic parameters, 
as displayed in Fig. [3] In the left panels Tq (a), T[u-b] (c) and 
T[e\] (e) are compared to T^ff and in the right panels log^g 
(b), log g[u-h] (d) and log^^ (f) are compared to \ogg. Dotted 
lines indicate internal uncertainties like in Table |2] for calibra- 
tions independent of the luminosity class (Alog gQ- 0.08 dex, 
Alogg[u-h]- 0.19dex and AT[c\]- 1060K). Otherwise, they 
indicate averaged Icr values for different luminosity classes 
weighted by star number (ATq= 700 K, AT^^-b]^ 1040K and 
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Table 2. Internal and external uncertainties for temperature and 
gravity determinations using the calibrations from Sects. 1471114. 31 



ibrations, as there are more stars with measured UVB indices 
than with other filters. 





AT 


Int. 


Ext. 


Alogg 


Int. 


Ext. 


V 

IIIJV 


ATq: 
ATq: 


738 K 
385 K 


800 K 
530 K 


Aloggg: 
Aloggg: 


0.08 
0.08 


0.23 
0.17 


IV,V 

III 




1060 K 
960 K 


300 K 
280 K 


Alogg[„-6j: 
AIogg[„_6]: 


0.19 
0.19 


0.02 
0.04 


V 

mjv 




1060 K 
1060 K 


290 K 
290 K 


Alogg/,: 
Alogg/,: 


0.07 
0.08 


0.01 
0.01 



A log gj3- 0.08 dex). Dashed lines indicate perfect agreement and 
symbols are encoded as in the legend. Stars #4 and 26 were ex- 
cluded from the fits in Fig. |3f as they are too evolved (i.e., be- 
yond core H-burning) in comparison to the rest of the sample. 

As indicated in Table|2] the internal uncertainties based on Q 
are smaller than in the other cases, but the external uncertainties 
are larger In further comparisons to other work, we take as ref- 
erence values the internal uncertainties based on Q: for tempera- 
ture ATq= 700 K (averaged value for different luminosity classes 
weighted by star number) and for gravity A log gQ- 0.08 dex, as 
in Fig. [3] This is because our special interest in the 2-based cal- 



5. Our spectroscopic parameters vs. other works 

5.1. Photometric Teff calibrations 

Our Tg-scales derived from multiple ionization equihbria are 
compa red in Fig. |4^ to photometric calibrations of Daflon et alj 
(Il999h for luminosity class V and of lLvubimkov et al . (2002) for 
\ogg =3.48 and 4. 02. In add ition, our r[„_i]-scales are compared 
to that of Napiwot zki et alJ (il993) in Fig.|4j5. Our scales tend to 
be hotter at higher temperatures than the photometric scales. 

Our Tq for luminosity class V (upper solid curve, blue 
in the online v e rsion) is ~2000K hotter than that proposed 
by iDaflon et alj (|1999) at the highest temperature, the differ- 
ence decreasing until both scales cross each other at ~21 000 K, 
and remain simi l ar at shghtly lower T^ff. The comparison to 
Lvubim kov et al.l (|2002|) has to be done more cautiously because 
our sample has a variety of surface gravity values. Only four 
stars have similar gravities to 4.00: #3 (III), 7 (IV), 16 (IV) and 
29 (V). We compare their upper curve (dashed dotted line) to our 
lower curve (solid line, red in the online version). Our relation is 
~2000 K hotter than theirs at the highest temperatures and both 
practically coincide below 2 1 000 K. 

Figure |4j3 sho ws the comparison of ou r T[u-h] to the calibra- 
tion proposed bv iNapiwotzki et aP (Il993h in their Eqn. 9. Their 
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and our curve for luminosity class III practically coincide below 
19000K, and the smallest discrepancy to our IV, V relation is 
~1500K at Tetf ~18 000K. In contrast, our scales are steeper 
than their curve, resulting in larger differences at higher Teff, 
with maximum values of ~1500K for luminosity class III and 
~6000K for IV, V. Possible reasons for these discrepancies are 
discussed at the end of next section. 



5.2. Other determinations for individual stars 

Besides the photometric temperature scales discussed in 
Sect. 15.11 other determinations of temperature and surface grav- 
ity are found in the literature for several stars of the sam- 
ple. Whenever possible, we also compare values of microtur- 
bulence as all parameters are interrelated. Figures |5]|9] illustrate 
the differences between the atmospheric parameters derived in 
other works and our spectroscopic Tgff (a), \ogg (b) and ^ (c). 
Dashed lines indicate full agreement between other determina- 
tions and our values, and dotted lines in Figs.|5]|9^,b denote ATq 
(700 K) and Alog^Q (0.08 dex), like in Fig. [3] as a reference 
to our photometric scale. Note that uncertainties in our spectro- 
scopic parameters are even smaller than the Icr-errors for the 
g-calibrations: ~300K for T^ff and ~0.05 dex for \ogg. 

It is beyond the scope of this paper a detailed descrip- 
tion of other analysis methods and the causes of discrepancies 
between their and our r e sults. Based on our previous studies 
(iNieva & Przvbillal l2007l l2008l 1201 Oa'U. we propose possible 
reasons that might cause systematic offsets. However, we should 
bear in mind that in most cases the discrepancies are likely due to 
a combination of various effects, which are difficult to be quan- 
tified beforehand. 

iKilian et^ jl99l|) . This work presents a spectroscopic deter- 
mination of gravity and effective temperature of stars in the so- 
lar neighborhood from NLTE hydrogen line profiles (Hj6, Hy, 
H6) and the silicon ionization equilibria Siii/iii and/or Siiii/iv 
(depending on the stellar temperature). Their procedure is simi- 
lar to ours, based on an earlier version of the codes Detail and 
Surface, however with fewer spectroscopic indicators than those 
shown in Table IBT] Their silicon abundances are based on the 
model atom by Becker & Butler (1990fl like in our case, but 
on model atmospheres by Gold ( 1984) with much less line blan- 
keting than ATLAS9. Their microturbulences are derived in the 
standard way from On lines. Figure|5]displays comparisons be- 
tween their and our stellar parameters for 1 1 stars in common, 
identified like in Table [1] A clear trend is noticed in Fig. 
for effective temperatures, being their values ~2000K higher at 
Teff -20 000 K and -2000 K lower at T^s -34 000 K. Less blan- 
keted model atmospheres affect the line-blocking and hence the 
temperature structure of the stellar atmosphere. This could be a 
reason for discrepancies with our temperatures. Furthermore, for 
cooler stars, the use of Si n in the ionization equilibrium can also 
cause an offset in the derived temperature. In contrast, their grav- 
ities agree with our spectroscopic values and lie within the Icr 
uncertainty of loggg for most stars (Fig.|5]3). Their microturbu- 
lences tend to be larger than ours (Fig.|5};), reaching a maximum 
difference of ~6kms 



We do not analyse Si ii lines because tlieir syntlietic spectrum pro- 
vides erroneous line intensi ties, as reported in lSimon- Diazl ilOVi) and 
iNieva & Simon- Diazl ( l201lh . 



ICunha & Lamberii jl992|) . The parameters determined in this 
work were used in a series of subsequent studies of stars in the 
Orion region bypu nha & Lambert ( 1994), Cunha et al. (2006|), 
iLanz et al.1 (l2008ir and IDaflon et all fl2009). Their approach for 
deriving effective temperatures and surface gravities is based on 
an iterative scheme usi ng older calibr a tions of Stromgren pho- 
tometry modified by Gies & Lambert (1992) coupled with fits 
to the wings of Hy profiles computed in LTE (Kurucz 1979'). 
Microturbulences were determined from On lines in NLTE. 
Figure|6]displays the differences between their and our parame- 
ters for 1 3 stars in com mon. In spite this sa mple differs from that 
in the comparison with Kilian et al.l (1199 ll) , we notice in Fig.|6^ 
an overall similar trend in temperature than in Fig |5h, except 
for star #22. But in this case, there are more stars with tempera- 
tures within our Icr uncertainty margins for Tq. It is difficult to 
identify the sources of discrepancies in temperature. They could 
be related to the details of the photometric calibrations, i.e. the 
use of less blanketed model atmospheres to compute the syn- 
thetic photometric indices and to the fits to Hy in LTE in the 
iterative procedure. Concerning the gravities (Fig.|6}3), their val- 
ues are higher (up to 0.4 dex) than our spectroscopic gravities 
for most stars, lying above the Icr uncertainty of our \og gg cal- 
ibration. An exception is the star #21, that shows a difference 
of -0.4 dex to our value. The discrepancies in surface gravities 
could be related to the use of LTE synthetic profiles of Hy, that 
have been shown to give larger gr avities than those computed 
in NLTE dNieva & Przvbillall2007h . Since their approach is it- 
erative, this also would affect the determination of temperature. 
However, differences larger than +0.2 dex cannot be explained 
with NLTE effects on the Balmer lines wings. The continuum 
normalization of Hy could also play a role in this respect, since 
this is an important sour ce of uncertai nty in surface gravity deter- 
minations (Ni eva & Przvbillall2010allbh . Their microturbulences 
are larger - by 5 to 10km s~' - than those derived b y us. 

In contra st, in Cunha et all (l2006h . iLanz et al.1 (2008) and 
IDaflon et all ([2009 f| the stars # 1 1, 21, 22, 27 and 28 were ex- 
cluded from the studies, which lead to maximum differences to 
our work in Teff of up to lOOOK and in \ogg of -1-0.05-0.3 dex. 
I.e., the selected sub-sample of stars analyzed in those works 
presents a better agree ment with our results t han the original 
complete sample from Cunha & Lambert! (Il992h . however a sys- 
tematic offset in gravities and microturbulences still persists. 

iMorel & ButleH J2008V This work presents an interesting com- 
parison between atmospheric parameters derived in a standard 
iterative spectroscopic way, i.e. via NLTE Si ii/iii and/or Si iii/iv 
ionization equilibria (Morel et al. 2006) and an new Nei/ii ion- 
ization equilibrium. In addition, the fits of NLTE profiles to 
He, H6, Hy and Hfi are used for surface gravity determina- 
tion. The computations are performed with similar versions of 
Detail and Surface and the same model atmospheres than those 
used by us. For Si, a model at om by D. J. Lennon is em- 
ployed dTrundle & Lennoiil l2005h while for Ne a new model 
by K. Butler is applied for the first time. Temperatures derived 
with Ne are cooler by 825 K than those derived with Si. In 
both cases, microturbulences were determined from On lines. 
It is worth noticing that they have analyzed 3 stars that were 
excluded from our first sample (Paper I) because they turned 
out to be double-lined spectroscopic binaries, namely HD 35468 

^ Pa rameters for 2 stars in common iiave been revised in lDaflon et al.l 
l l2009h via Sn/m ionization equilibrium. For star #1, offsets of +530 K 
and +0.05 dex and for #24, -480 K and -0.05 dex in temperature and 
surface gravity, respectively, have been applied. 
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Fig. 7. Like Fig. |5] but for parameters derived bv iMorel & ButleJ (l2008h via ionization equilibria of Ne i/ii (diamonds) and Si n/m- 
Si iii/iv (triangles). 
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Fig. 9. P anels (a) a nd (b) like Fig.|5j but for parameters derived with the calibrations of' Moon & Dworetskvl ( |1984|) . later improved by 
iNapiwotzki et alJ ( 11993) . Panel (c) shows the difference between the [u - b]-hased temperature scale proposed bv lNapiwotzki et aP 
(Il993h and our spectroscopic values. Diamonds indicate stars with Teff < 30 000 K and squares with T^ff > 30 000 K. 



(HR1790, Bellatrix), HD 163472 (HR6684, V2052Oph) and 
HD 214993 (HR 8640, 12 Lac). Our samples share 6 stars, which 
parameters are compared in Fig. I?] Overall, we notice a trend 
in temperature (Fig. |7^) in a similar way then the two previous 
cases, the maximum differences being ~1000 at ~21 OOOK and 
-2000 K at -27 000 K. 

It is difficult to evaluate the causes of discrepancies in the 
results. It is not clear whether the new Si ii mo del atom provides 
impr oved line intensities than the model by iBecker & Butled 
(1199 0). If the problem persists, it would affect temperatures 
lower than ~23 000 K (the exact limit depends o n the S/N of the 
spectr um). In addition, the Ne model applied in lMorel & Butler! 
(l2008h showed difficulties to achieve Nei/ii ionization equi- 
librium and to match the temperatures derived with Siii/iii/iv 
ionization equilibria. This Ne model atom wa s subsequently 
impro ved with more accurate log gf values in iPrzvbilla et alj 
(1200 8") resulting in consistent results between Ne i/ii and other 
ionization equilibria: Hei/ii, Cii/iii/iv, Oi/ii, Siiii/iv, removing 
discrepancies found previously; see also Table lB.ll of the present 
paper. Furthermore, their surface gravities (Fig. |7}5) are ~0.1- 
0.45 dex lower than those found by us. These offsets are unex- 
pected, as their analysis is based on the same codes and a sim- 
ilar H model than those employed in our work. An offset in the 
temperature determination is not sufficient to explain such large 
differences in log g. The remaining sources of systematics to be 
investigated in tho se results is the continuu m normalization of 
the Balmer lines (Nieva & Przvbillall2010allbl) . that might cause 
differences larger than 0.2 dex in the derived surface gravities. 
Finally, there is a hint of a trend in microturbulences, except for 
star #6 (Fig. E};). 

ISimon-DfazI 1201 Ol) . Stellar parameters were determined in 
this work spectroscopically by fitting the Balmer lines and 
matching He i/ii and/or Si iii/iv ionization equilibria. The spec- 



trum synthesis was computed with the unified atmosphere code 
Fastwind (Puis et al. 2005), that models the stellar atmosphere 
and wind simultaneously in NLTE. The same observed spec- 
tra than those analyzed in Paper II have been analyzed in 
ISimon-D faz (2010). Despite the different atmospheric structure, 
input atomic data and line-formation computations, the temper- 
atures and gravities derived in this work agree overall very well 
with our determinations, as can be seen from Figs.[8^,b, i.e., the 
differences lie within the Icr-uncertainty of the g-calibrations 
for most stars. Microturbulences also show a good agreement, 
considering the different computations, reaching a maximum 
difference of ~4 km s ' (Fig. [8]:). 

INapiwotzki et al.l j1993|) . In Sect. [H our Trtz-fci calibratiori s 
were compared with that proposed by Napiwot zki et aP (Il993h . 
In Fig. |9^,b, our spectroscopic parameters are compared to 
tem peratures and gravities c omputed with an improved version 
by Napi wotzki et al.l (Il993h of a code by iMoon & Dworetskvl 
(1984). Fig. |9]; shows a comparison of our spectroscopic tem- 
peratures to those derived with the [u - /?] -temperature scale 
bv Napiwotzki et aP (Il993h . The code uses as input Stromgren 
ubvyP indices and interpolates within a grid of synthetic colors 
iteratively until convergence in effective temperature a nd sur- 
face gravity is reached. The modification by Napiwotz ki et al.l 
(Il993h on the gravities is based on fits to Balmer lines. The 
code provides a final temperature and gravity, hereafter Tai93 and 
loggN93, respectively. It also returns a temperature from their 
calibration to the [u - /?] -parameter, hereafter TNgsiu-b], that cor- 
responds to their Eqn. 9, discussed in Sect. l5.1l (see also Fig.|33). 

Figure|9^ displays the difference between their rN93 and our 
Teff. There is a very good agreement for lower temperatures, up 
to T^ff- 24 OOOK. At higher temperatures there is a larger scat- 
ter, their values being smaller than ours from a few hundred up 
to ~3000K. Exceptions are star #12, with perfect agreement. 
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and star #22, ~2500 K hotter in their approach. We recall that a 
non-standard reddening law is found for this star, as it is located 
in a dense nebulous region, therefore a photometric determina- 
tion of its parameters is dubious. The vertical dashed-dotted line 
in Fig. corresponds to the upper temperature limit to which 
they expect their determination to be valid. Thick squares denote 
stars with reff> 30 000 K and diamonds the rest of the sample; 
this corresponds to the upper limit to which they expect their 
method to be valid. From Fig.|9^ we notice, however, that their 
upper limit is likely reff~ 24 000 K. Figure |9}3 shows the dif- 
ferences between their gravities, loggN93, and our \ogg. There 
is a good agreement between both methods for several stars. A 
few objects also show larger discrepancies up to +0.3 dex for 
stars denoted with diamonds and, as expected, even larger for 
those marked with squares, with temperatures beyond their up- 
per limit (+0.6 dex for star #30 and -0.5 dex for star #21, not 
shown in Fig.|9j3 because they lie outside the plot range). Their 
rN93[u-bj, extracted from the same code, is compared with our 
spectroscopic temperature in Fig.|9};- The differences to our work 
show a similar behaviour than Fig.|9^, but with a larger scatter. 
Stars #2, 21 and 30 are not shown because their differences are 
larger than -3000 K. Following the comparisons in Fig. |9l we 
recommend to adopt Tnqj for Teff < 24 000 K when using their 
approach. 

Several reasons could lead to discrepancies between their 
and our results at higher temperatures. We note that only four 
stars hotter than 24 000 K are present in their calibration sam- 
ple to derive the rN93[u-b] -scale. Three of these stars have lu- 
minosity class III, therefore the rN93[u-b]-scale underestimates 
temperature of stars with luminosity class V (most objects in 
this work). Two of these four stars are present in our sample: #3 
and #19. rN93[u-b] for star- #3 (HD 63922, HR3055) is 1200 K 
lower than our value, and that for star #19 (HD 34816, HR 1756) 
is 2800 K than our temperat ure. Furthermore , their calibration is 
based on temperatures from Malagni ni et al.l (11986), determined 
from comparisons of obser ved spec t ral ene rgy distributions to 
fluxes computed in LTE bv lKuruc3 (Il979l) . Hence, discrepan- 
cies with our results could also be caused by less line blanket- 
ing in their model atmospheres, neglection of NLTE effects in 
the fluxes, and also in the computation of synthetic indices and 
B aimer lines. 



Table 3. Comparison of T^ff- and BC-scales from this work to 
literature data. 



Sp. T # Teff logg BC TsKi2 J^CsKSl Tqc94 Tcoq 



BlIV 14 27000 4.05 -2.67 

B 1.5 IV 16 23000 3.95 -2.32 

8 22000 3.85 -2.18 

B2IV 7 22000 3.95 -2.20 

18 21250 3.80 -2.11 

17 20750 3.80 -2.05 

B3IV 20 17500 3.80 -1.65 

18 21250 3.80 -2.11 

B0.2III 3 31200 3.95 -2.92 29000 -2.88 

B0.5III 12 30000 4.05 -2.88 

B3III 15 15800 3.75 -1.38 17100 -1.60 

SK82: 'Sc hmidt-Kale3 ( fl982h . GC94: iGrav & Corbaihl fT994h . COO: 
Cox ( 200dv 
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5.3. Our recommended parameter determination 

Our scales derived spectroscopically tend to be hotter 
at higher temperatures than the photometric scales pro- 
posed bv [Pafl on et all (Il999h . iLvubimkov et al.l (l2002h and 
iNapiwotzki et al.l (11993^ ). Furthermore, our values also tend to 
be higher at Teff > 25 000 K and lower at lower temp e rature s 
when compared with va lues derived by Kilianet al.1 (Il99lh . 
ICunha & Lambert! ( Il992h and lMorel & Butled ( l2008h " using dif- 
ferent analysis methods, synthetic spectra and obs ervations. In 
cont rast, we note that temperatures determined by ISimon-DiazI 
(12010.) using the same set of high-quality spectra than those anal- 
ysed by us, other state-of-the-art model atmosphere and spectral 
synthesis in NLTE and simultaneous silicon ionization equilibria 
and fits to Balmer and helium lines, render values in agreement 
to our determinatio ns. The excellent agreement between our and 
ISimon-Dfad (12010 ') temperatures lead us to conclude that our 
temperature determination based on high quality data, state-of- 
the-art spectral synthesis in NLTE and several independent tem- 
perature indicators are reliable and almost free of large remain- 
ing systematic effects. If adopting significantly different values 
in the higher temperature domain, like suggested by previous 



studies, the ionization equilibrium of no species will be matched 
any longer, hence the global and detailed observed spectra will 
not be reproduced by the synthetic model. 

We recommend, when possible, to perform a parameter de- 
termination like the full spectroscopic r nethod prese n ted by us 
in previous works (e.g. Paper I) or by ISimon- Dfa3 (120101). If 
the observed spectra or the spectral synthesis are not available, 
a good parameter estimation can be provided by our scales pre- 
se nted in Sect. [Hof this wor k. For reif< 24 000 K, the scale rN93 
by INapiwotzki et al.l (1 19931) also rend ers reliable te mperatures, 
and for T^s^JAQOOK, the scales bv lDaflon etaTI (11999 ) and 
Lvubimk ov et al.l (l2002l) also give reliable results (note the lower 
limit indicated by each scale). 



6. Spectral type calibration 

In App.lAl we re-evaluate the spectral classification of the sam- 
ple stars. Slight revision of spectral sub-types and luminosity 
classes have been performed in several stars, as indicated in 
Table [T] The relationship between spectral types and luminos- 
ity classes and the stellar parameters Tetf , log g as well as their 
corresponding value of bolometric correction (BC) is summa- 
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4.6 4.5 4.4 A. 6 4.2 4.1 4.6 4.5 4.4 4.5 4.2 4.1 

log (Teff) log (Teff) 

Fig. 10. Left panel: bolometric corrections vs. effective temperature as derived for the sample stars. A typical error bar is shown, the 
line indicates a regression curve to the data. Right panel: comparison of our T^s-BC relation with data from the literature. See the 
text for details. 



rized in Table[3] There, also recommended values of temperature 
and /or B C from Schmidt-Kaler ( 1982), Gray & Corbally ( 199^ 
and lCoxl fcOOO ) are given for comparison. 

In cases where more than one star belong to a certain spec- 
tral/luminosity class, instead of adopting an average value for 
the parameters, we prefer to indicate the exact parameters of the 
stars to see how much can they variate. Stars #4 and 26 have 
been excluded from the calibration because they seem to be too 
evolved in comparison to the rest of the objects, as explained 
above. Our star sample is not large enough to cover all possi- 
ble combinations of spectral and luminosity classes within the 
range 09-B3 V-lIl, however the carefully derived data allow us 
to anchor a spectral type calibration to empirical values of tem- 
perature for luminosity class V and partially also for luminosity 
classes 111 and IV. Our calibration tends to be hotter than pre- 
vious reference values in most cases. This is in accordance with 
our photometric temperature calibrations, that are in general hot- 
ter than previous work. This table is also very useful to estimate 
stellar parameters once the spectral and luminosity classification 
are well known. 



7. Bolometric corrections vs. temperature 

Bolometric corrections play an important role in the determi- 
nation of luminosities and indirectly also of other fundamental 
stellar parameters computed from atmospheric parameters Teff 
and \ogg, e.g. radii, masses. Fundamental parameters can be di- 
rectly compared to predictions from evolution models in order 
to learn more about the physics of the stars. 

Taking advantage of our accurate atmospheric parameters 
determined spectroscopically for the present sample, we derive a 
semiempiric relationship for BC as a function of T^g. Bolometric 
corrections are computed from spectral fluxes calculated with 
ATLAS^fl for each star of the sample using their corresponding 
values of Teff and log g. As expected, a clear relationship is found 
for BC vs. Teff (Fig.[TOl left panel), practically independent from 
the luminosity class or the surface gravity. Figure[TOl right panel, 
illustrates a comp arison of our relationship to other work by 
IVacca etal.1 (ll996h .lCox (2000), Hubeny & Lanz (2003, 2007), 



Atmospheric structures and fluxes of OB main sequence and giant 
stars computed in LTE a re e quivalent to the full N LTE approach. See 
iNieva & Przvbillj j2007h and lPrzvbilla et all j201lh for details. 



and iMartins e t al.l (l2005 l) derived mainly from synthetic mod- 
els. O ur slope is less steep than those by Vacca et al. ( 1996) and 
IMartins et al] ( l2005h . our calibration extending to lower temper- 
atures . In add ition, our fit agrees with BC from iHubenv & Land 
(I2003L 12007 1) and for temperature s higher than -32 000 K also 
with those by Martins et al. ("2005). In contrast, reference values 
for BC recommended by Cox (2000) are systematically lower 
than our data. 

From the relationship found in Fig. [10] (left panel), we pro- 
pose a new calibration of the bolometric correction as a function 

of Teff: 



BC = 21.00-5.34 log(reif); Itr = 0.01 (10) 

The boundaries are limited by the observational data: 
15 800 < Teff < 34 000 K. For higher te mperatures, we recom- 
mend to use bolometric corrections from lHubenv & Land (l2003h 
or lMartins et alJ (l2005h . 

8. Summary and conclusions 

Based on comprehensive previous studies which lead us to deter- 
mine effective temperatures and surface gravities of 30 OB-type 
(09-B3 V-III) stars at high accuracy and precision, here we in- 
vestigate whether such parameters can be used to propose new 
calibrations to some photometric indices. The spectroscopic pa- 
rameters were derived via multiple ionization equilibria and si- 
multaneous fits to the Balmer lines based on our most recent 
spectrum synthesis for NLTE computations and high-quality ob- 
served spectra. We propose very useful relationships between 
atmospheric parameters and various photometric indices from 
the Johnson and the Stromgren systems. Effective temperatures 
can be estimated at a precision of ~400 K for luminosity classes 
III/IV and ~800 K for luminosity class V. And surface gravities 
can reach internal uncertainties as low as ~0.08 dex when using 
our calibration to the Johnson Q-parameter. Similar precision is 
achieved for gravities derived from the /3-index and the preci- 
sion decreases for both atmospheric parameters when using the 
Stromgren indices [cl] and [u - b\. In contrast, external uncer- 
tainties are larger for the Johnson than for the Stromgren cali- 
brations. Our uncertainties are smaller than typical differences 
among other methods in the literature, reaching values up to 
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+2000 K for temperature and +0.25 dex for gravity, and in ex- 
treme cases, +6000 K and ±0.4 dex, respectively. We had further 
re-examined the luminosity and the spectral classes for the star 
sample and also propose a parameter calibration for sub-spectral 
types. Furthermore, we present a new bolometric correction re- 
lation to temperature based on our empirical data, rather than on 
synthetic grids. 

The photometric calibrations presented here are useful tools 
to estimate effective temperatures and surface gravities of non- 
supergiant OB stars in a fast manner. Caution has to be taken for 
undetected double-lined spectroscopic binaries and single ob- 
jects with anomalous reddening-law, dubious photometric quan- 
tities and/or luminosity classes, for which the external uncertain- 
ties may increase significantly. We recommend to use these cal- 
ibrations only as a first step for parameter estimation, with sub- 
sequent refinements based on spectroscopy. 
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Appendix A: Spectral types and luminosity classes 

The photometric calibrations discussed in Sect. |4]depend on lu- 
minosity classifications of the star sample. Here, we re-examine 
their luminosity and spectral classification. As reference, we 
adopt the classification of one anchor point of the Morgan- 
Keenan system (star #11) and five primary standards from N. 
Walbo rn (#6, 7, 11, 19 and 23), all Hsted in Gray & Corbali^ 
(l2009l) . In a second step, spectral types and luminosity classes of 
the sample stars are examined via inspection of Balmer line s and 
sihcon/heUum ra t ios, as recommended by Wal bornI (|1971|) and 
iGrav & Corballvl(l2009h . Spectral type classification for the hot- 
ter stars (O9-B0.7) rely on the ratio Si iv^4089/Si iii/14552 and 
for the cooler stars (B1-B3) on Siii/i4128/Siiii/14552. In addi- 
tion, luminosity classes are tested via Si iv/l41 16/He 1/14121 for 
the hotter and Si m/14552/He 1/14387 for the cooler stars. 

Figure lATI displavs the relation of line ratios as spectral type 
and luminosity class indicators for the hotter (upper panel) and 
the cooler stars (lower panel). Only few cases - stars #2, 8, 10, 19 
and 26 - required a coiTection in luminosity class. Furthermore, 
slight coiTections to (sub-)spectral types for objects #2, 5, 8, 10, 
13, 16, 28, 29, have also been implemented, as listed in Table[T] 
The revisions are facilitated by the higher quality of our spectra 
in comparison to older photographic spectrograms used in the 
past, affecting in particular the weaker silicon lines. Star #15 is 
too cool, i.e. has a very low Si ii/l4128/Si iii/l4552 ratio, therefore 
is excluded from Figure |AT| (lower panel). 

Appendix B: Spectroscopic indicators for 
temperature and surface determination 

Table IB. II indicates which ionization equilibria in combination 
with fits to several hydrogen lines, have been matched in order 
to determine effective temperatures and surface gravities for the 
sample stars. Dots denote fits to spectral lines and boxes indicate 
cases where ionization equilibria were matched. A row denotes 
the simultaneous fits to different elements/ionization stages us- 
ing the same set of atmospheric parameters. 
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M. F. Nieva: Temperature, gravity and bolometric correction scales 
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Fig. A.l. Ratios of equivalent widths for temperature (abscissa) 
and gravity-sensitive (ordinate) line pairs as basis for the spectral 
type and luminosity classification of the star sample (Table [TJ. 
Upper panel: hotter stars O9-B0.7; lower panel: cooler stars BI- 
BS. 



Table B.l. Spectroscopic indicators for T^g and logg determi- 
nation. Dots indicate spectral fitting and boxes denote ionization 
equilibria. 

HD Teff H HeiHeii CiiCiiiCiv OiOii NeiNeii SiiiiSiiv FeiiFeiii 





1 


36591 


27.0 


2 


61068 


26.3 


3 


63922 


31.2 


4 


74575 


22.9 


5 


122980 


20.8 


6 


149438 


32.0 


7 


886 


22.0 


8 


29248 


22.0 


9 


35299 


23.5 


10 


35708 


20.7 


11 


36512 


33.4 


12 


36822 


30.0 


13 


36960 


29.0 


14 


205021 


27.0 


15 


209008 


15.8 


16 


216916 


23.0 


17 


3360 


20.7 


18 


16582 


21.2 


19 


34816 


30.4 


20 


160762 


17.5 


21 


37020 


30.7 


22 


37042 


29.3 


23 


36959 


26.1 


24 


37744 


24.0 


25 


36285 


21.7 


26 


35039 


19.6 


27 


36629 


20.3 


28 


36430 


19.3 


29 


35912 


19.0 


30 


46202 


34.1 
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